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I  INTRODUCTION 


The  Institute  for  Telecommunication  Sciences — Naval  Research 
Laboratory  (ITS-NRL)  area-coverage  model  has  been  developed  to  serve 

three  purposes— first,  the  analysis  of  measured  area  coverage  and  radar 
backseat ter  data;  second,  the  prediction  of  the  performance  of  communi¬ 
cation  and  radar  systems;  and  third,  as  a  tool  to  be  used  in  frequency 
and  scan  management  of  OTH  radar  systems.  The  development  and  use  of 
this  model  for  these  purposes  have  been  discussed  in  many  previous 
reports.1”3 

The  purpose  of  this  paper  is  to  describe  the  present  status  of 
the  model  and  to  illustrate  its  use  in  analyzing  radar  backscatter  data. 


II  MODEL  DESCRIPTION 

The  model  is  divided  into. two  parts;  the  first  half  provides 
a  description  of  the  ionosphere  along  *my  look  direction  of  a  radar, 
given  its  location  (Figure  1).  If  ionospheric  data  are  already  available 
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either  from  previous  analysis  or  measurement ,■ the  first  half  of  the  pro¬ 
gram  may  be  skipped  and  the  data  entered  directly  into  the  second  half 
of  the  model  (see  below  for  details)  If  no  data  are  available,  the  ITS 
world  maps  of  ionospheric  indices  are  used  to  generate  the  description  of 
the  ionosphere.4-9  These  maps  have  been  in  use  several  years  except  for 
the  Fl-layer,  which  is  described  in  detail  in  Ref.  9  and  summarized  in 
the  Appendix.  The  ionosphere  is  described  initially  by  four  ionized 
regions — the  regular  E,  the  FI  if  present,  the  F2,  and  the  sporadic  E. 

The  E,  FI,  and  F2  are  described  each  by  four  indices — the 
critical  frequency,  the  height  of  maximum  ionization,  the  semi-thickness 
of  the  layer,  and  an  absorption  factor  (described  in  detail  below,  it 
corresponds  roughly  to  the  averaged  effect  of  the  product  of  collision 
frequency  and  electron  density).  The  sporadic  E  is  described  by  four 
parameters — the  height  of  reflection,  the  lower  decile,  the  median,  and 
the  upper  decile  value  of  median  critical  frequency.  These  sixteen 
indices  form  the  first  rough  description  of  the  ionosphere;  if  desired 
these  may  be  entered  as  input  data  rather  than  generated  from  the  world 
maps. 

The  second  step  of  the  model  is  to  generate  an  electron-density 
profile  describing  the  E,  FI,  and  F2  portions  of  the  ionosphere.  At  each 
selected  height  the  square  of  the  plasma  frequency  is  generated  using  the 
equations  for  a  parabolic  layer  (Figure  2): 

4  =  ^ci-(~s)2]  a) 
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(approximately  110  km)  is  used  to  start  the  electron-density  profile  at 
the  bottom.  This  yields  virtual-height  ionograms  and  raypaths  that  agree 
well  with  measured  data.  If  desired,  an  option  exists  to  tail  off  the 
E  layer  with  a  linear  layer  down  to  a  desired  level.  For  example,  if 
the  bottom  of  the  E  layer  is  taken  to  be  90  km,  a  linear  tail  to  the 
E  layer  starting  at  100  km  and  extending  down  to  70  km  will  produce  a 
virtual  height  of  110  km  at  a  true  height  of  110  km  [f^  =  U.75(fQE)2 

for  a  semithickness  of  20  km].  Again,  it  is  optional  to  enter  any 
electron-density  profile  desired.10 

The  fourth  step  of  the  model  is  to  generate  a  virtual-height 
ionogram  from  the  electron-density  profile  (Figure  3) .  Here  a  numerical 
integration  is  used  to  evaluate  the  virtual  height,  h/: 


hr 

h7  =  J  ji'OOdH 


where 

I j/  =  Group  index  of  refraction 
y  =  Ratio  of  plasma  frequency  to  critical  frequency 


h 


2 

Height  corresponding  to  f^  , 


In  the  program  a  table  of  y/  as  a  function  of  (Figure  3, 
lower  right)  is  stored.  For  each  value  of  n^,  the  graph  of  y/  (Figure  3, 
upper  left)  is  generated  using  the  table  of  | y/  and  y^.  The  virtual 
height  corresponding  to  hr  is  then  the  area  under  the  y/vs  h  curve  (Fig¬ 
ure  3,  upper  left).  This  procedure  gives  a  virtual-height  ionogram 
(Figure  4)  with  a  minimum  of  calculations  for  any  tabulated  electron- 
density  profile.  The  virtual-height  trace  completes  the  description  of 
the  ionosphere  except  for  the  absorption  factor  discussed  in  detail  be¬ 
low.  This  completes  the  first  half  of  the  program  which  can,  of  course, 
no  used  as  a  complete  program  to  generate  ionospheric  indices,,  electron- 
density  profiles,  and  virtual-height  ionograms  on  a  worldwide  basis  from 
the  ITS  median  prediction  coefficients. 

The  second  half  of  the  program  is  designed  to  specify  the  area 
coverage  of  a  transmitter  operating  in  the  decameter  wave  band.  The 
expected  input  is  virtual-height  ionograms  describing  the  ionosphere  at 
selected  sample  points  along  the  desired  azimuth.  This  input  can  either 
be  that  generated  by  the  first  half  of  the  model,  or  it  can  be  entered 
separately  from  data  obtained  elsewhere.  The  following  simplifying 
assumptions  will  be  made  to  derive  the  raypath: 

(1)  Horizontal  and  azimuthal  variations  in  the  ionospheric 
electron-density  profiles  will  be  ignored  (for  each  hop;  on 
a  multihop  path  different  sample  areas  are  used). 

(2)  Magnetic-field  variations  are  ignore 

(3)  The  ionosphere  is  spherically  symmetrical. 

With  these  assumptions  the  raypath  is  completely  specified  by  Bouger's 
rule,  which  is  equivalent  to  Snell's  law  for  raypaths  v/ithin  a  spheri¬ 
cally  symmetric  medium  (Figure  5),  and  is  defined  by: 


(irsinQ  =  constant 


(3 


f,,  =  f  sec0  (3 

O  V 

Index  of  refraction 

Radtus  from  the  center  to  the  fay 

Anglo  between  the  ray  and  the  normal  to  the  sphere 

Angle  between  the  apparent  raypath  and  the  normal  to  the 
sphere  at  the  true  height  of  reflection 


/  i 


FIGURE  5  RAYPATH  GEOMETRY 


(4) 


sin0 


sinA 

1  +  Vre 


A  =  Take-off  angle  of  the  ray 

hp  =  True  height  of  reflection. 

r  =  Earth's  radius, 
e 

The  model  used  Eqs.  (4)  and  (3b)  to  generate  Table  1.  The 
first  column  is  the  table  index.  This  also  corresponds  to  ten  times  the 
vertical  sounding  frequency— e.g. ,  the  17th  row  of  the  table  is  for  a 
frequency  of  1.7  MHz.  The  second  column  is  the  true  height  of  reflection, 
hr;  the  third  column  is  the  virtual  height  of  reflection,  h';  and  the  col¬ 
umns  following  are  the  equivalent  oblique  operating  frequencies  of  the 
transmitter  for  ray  paths  with  corresponding  take-off  angles — e.g.,  the 
column  under  0  is  for  raypaths  with  a  take-off  angle  of  <f  .  Table  1  has 
a  one-degree  increment  between  take-off  angles  for  this  example;  the 
increment  is  arbitrary  for  the  program.  When  the  information  contained 
in  Table  1  is  plotted  in  constant-frequency  contours  of  virtual  reflec¬ 
tion  height  and  distance  as  in  Figure  6,  the  displayed  contours  are 
called  ref lectrices.11 * 13  At  each  given  operating  frequency  the  area 
coverage  is  found  by  interpolating  in  Table  1  for  the  desired  raypaths. 
This  procedure  yields  the  desired  reflectrix  for  the  given  frequency  in 
the  form  of  a  table  of  raysets  (the  firslf  six  columns  of  Table  2),  No 
calculations  are  necessary  at  each  frequency,  since  all  the  desired  in¬ 
formation  is  contained  in  Table  1.  For  two-hop  modes,  a  table  of  re¬ 
flectrix  information  (as  in  Table  1)  is  generated  for  the  second-hop 
sample  area,  and  at  the  given  operating  frequency  a  table  of  rayset 
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FIGURE  6  REFLECTRICES 


information  is  generated  (as  in  Table  2) ,  The  two-hop  modes  are  found 
by  matching  the  take-off  angle  of  the  second  hop  (DELI  as  in  Table  2)  with 
the  arrival  angle  of- the  first  hop  (DEL2  in  Table  2).  Note  that  there 
has  been  no  mention  of  individual  layers  (E,  FI,  or  F2)  since  the  electron- 
density  profile  was  generated.  The  ionosphere  is  treated  as  a  single 
region  by  the  program  and  all  possible  mode  combinations  are  generated.  In 
order  to  keep  the  traditional  layer  nomenclature,  the  raypaths  are  named 
according  to  where  their  equivalent  vertical  frequencies  lie  (e.g. ,  below 
f0E);  thus  the  modes  may  be  E-E,  F2-F2,  or  E-F2,  according  to  the  label 
of  each  hop.  Since  the  sporadic-E  modes  are  assumed  to  exist  with  some 
degree  of  probability  with  reflection  at  a  constant  height,  the  ray set  in¬ 
formation  is  the  same  for  each  operating  frequency  and  there  is  no  need 
to  generate  the  equivalent  of  Table  1  for  the  Eg  modes. 

The  model  then  uses  the  rayset  information  (Table  2)  to  add 
losses  to  each  raypath.  For  a  backscatter  radar,  the  radar  range  equution 
is  used  to  generate  the  si?nal-to-noise  ratio  for  each  raypath.  This  in¬ 
formation  is  then  used  to  generate  the  area-coverage  plot  in  Figure  7. 

The  noise  information  includes  atmospheric,  galactic,  and  man-made  sources. 
An  additional  man-made  noise  level  may  be  entered  as  part  of  the  input. 
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The  losses  are  those  as  calculated  by  Lucas  and  Haydon6  with  the  addition 
of  sporadic-E  losses  and  the  deviative  and  nondeviative  loss  factor,  The 
sporadic-E  losses  are  reflection  losses  for  modes  reflected  from  the  ER 
layer  and  obscuration  losses  for  raypaths  penetrating  the  Eg  layer.  The 
losses  depend  on  the  probability  of  a  sporadic-E  mode  existing  at  the 
given  operating  frequency  and  at  the  given  take-off  angle,  and  is  defined 
by 


P  =  Nor( 


f  -  fs, 


(5) 


where 


Nor  indicates  the  normal  distribution 


f  =  Operating  frequency 


f  =  Median  MUF  for  the  given  take-off  angle 


a  =  Standard  deviation  of  f  E  , 

o  s 


The  reflection  loss  is  then 


Lr  =  8.91(P) 


-0.7 


(6) 


and  the  obscuration  loss  is 


Lq  =  10  log10(l  -  P*) 


(7) 


The  ionospheric  losses  are  based  on  the  absorption  losses  for  D-  and  E- 
region  absorption  of  Lucas  and  Haydon.6 
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The  absorption  loss  in  the  quasi-longitudinal  approximation 


(Q.L.)  is 


rL(fvxf  -f  n 

L(f  )  =  - 1 -  .  ...  sec0 

L  -  fe  J 


(8) 


where 

L(fQ)  =  Absorption  for  the  oblique  path  with  operating 
frequency  f 

L(fy)  =  Absorption  for  vertical  paths  with  oper,  i;ing 
frequency  fv 

f  =  Longitudinal  component  of  the  electron  gyrofrequency 
c 

0  =  Angle  of  incidence  of  the  ray  at  the  absorbing  region 
(approximately  100  km) . 


The  term  in  brackets  in  Eq.  (8)  is  that  fitted  by  Laitinen  and 
Haydon  for  both  vertical  and  oblique  F-layer  paths,  13  Since  the  data 
did  not  include  frequencies  below  3  MHz,  the  equation  was  modified  by 
Lucas  and  Haydon5  to 


LA  = 


_(f  + 


677.2 


tj 


1.98 


+  10.2 


sec0 


(9) 


An  additional  loss  factor  for  E-layer  modes  has  been  derived: 


LE 


1.359  -  19.841 


(10) 


where 

fg  =  Median  MUF  for  the  E  layer  at  the  given  elevation  angle 
fQ  —  Operating  frequency. 
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An  additional  loss  term  is  now  included  that  takes  into 
account  the  deviative  absorption  effects  of  the  higher  layers  (of  most 
importance  for  high-angle  rays).  The  theoretical  loss  is 

L  =  P  k  X  dz  (11) 

V  J 

ray 

path 

where 

k  =  Propagation  constant  in  free  space 

n  =  (j,  -  ix  =  Refractive  index 

/ 

Neglecting  magnetic-field  effects  this  is 

\  -  f  (»'  -  V  <12> 

where 

v  =  Collision  frequency 
c  =  Speed  of  light 
h'  =  Virtual  height 
hp  =  Phase  height. 

In  order  to  obtain  an  expression  that  is  easily  calculated  and 
behaves  as  Eq.  (10),  we  approximate  the  bracketed  term  in  Eq.  8  with 

[  ]  =  A(fv)(h'  -  hr)  (13) 
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where 


h;  =  Virtual  height 

h^  =  True  ref lest ion  height 

A(£v)  =  An  empirical  fit  to  the  absorption  using  Eq.  (13). 

The  first  column  of  Table  1  is  the  index  corresponding  to  ten 
times  the  vertical  frequency;  the  second  column  is  the  true  height  of 
reflection,  hrj  and  the  third  column  is  the  virtual  height,  h/.  The 
fourth  column  is  the  loss  term  to  be  added  to  the  losses  calculated  by 
Lucas  and  Haydon5  [Eq.  (9)] — e.g. ,  the  sum  of  Eqs.  (10)  and  (13).  In  the 
program  A(fy)  is  now  a  given  constant  value  for  the  E  region,  the  E-Fl 
cusp,  the  FI  region,  the  F1-F2  cusp,  and  the  F2  region.  In  addition,  it 
is  optional  to  modify  this  term  by  any  constant  multiplication  factor 
(and  thus  eliminate  it  entirely  if  so  desired).  This  representation  is, 
of  course,  an  interim  step  until  the  completion  of  a  world  map  of  A(fv), 

III  DATA  ANALYSIS 

The  model  as  described  above  can  be  used  with  the  median  iono¬ 
spheric-index  maps  as  a  tool  to  predict  the  long-term  performance  of  HF 
radar  and  communication  systems.  Thu  use  of  the  model  as  a  tool  to  ana¬ 
lyze  and  update  the  predicted  median  data  will  nuW  be  described.  Assumed 
available  at  the  beginning  cf  the  analysis  are  the  median  ionospheric 
predictions  for  the  month,  including  the  yirtual-height  ionegram  (Fig¬ 
ure  4) ,  the  expected  backscatter  amplitude  versus  range  (Figure  7) ,  and 
the  antenna  pattern.  For  this  system  the  anter.na  is  composed  of  two 
horizontal  rows  of  elements  stacked  two  high.  The  antenna  patterns  re¬ 
sulting  from  feeding  these  two  rows  inphase  and  antiphase  are  shown  in 
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Figure  8.  A  scattering  coefficient  of  ct°  =  ~17  dB  was  used  to  describe 
the  sea  scattering,  properties  of  the  return. 

The  model  parameters  that  are  varied  are  those  relating  to 
electron-density  distribution— the  critical  frequencies;  heights  of  maxi¬ 
mum  ionization  and  semi-thicknesses  of  the  E,  FI,  and  F2-layers;  and  the 
critical  frequency  and  height  of  reflection  of  the  Eg  layer.  The  entire 
electron-density  profile  on  the  entire  h'  ionograms  can  be  modified 


FIGURE  8  ANTENNA  TRANSMIT-RECEIVE  GAIN  PRODUCT  GIVEN  AS  A 
FUNCTION  OF  RADIATION  ANGLES 


directly,  but  it  is  easier  to  consider  only  the  smaller  number  of  indices 
in  a  first  attempt  at  describing  the  actual  ionosphere. 

A  one-frequency  backscatter  sounding  is  sufficient  to  sample 
all  useful  heights  in  the  ionosphere;  but  to  determine  the  ion-height 
profile,  one  must  be  able  to  identify  tho  ground  range(s),  R,  or  radia¬ 
tion  angles  that  correspond  to  each  virtual  range,  Rf .  Positionable 
radiation- angle  directivity  would  provide  an  excellent  method  for  cou¬ 
pling  R  and  R*.  However,  in  this  exercise  the  tools  available  are 
limited  to  two  multi-lobed  antenna-radiation-angle  patterns  as  shown  in 
Figuie  8  and  a  chart  for  the  illuminated  region  as  shown  in  Figure  9. 

When  the  electron  density  versus  height  is  a  smooth  monotonic  function, 
the  antenna  lobe  structure  such  as  shown  in  Figure  8  enables  one  to 
easily  identify  R'  values  with  radiation  angles  and  then  compute  R  values 
and  heights.  Hence  a  more  complex  ionosphere  has  been  chosen  for  illus¬ 
tration. 

The  backscatter  level  plotted  against  time  delay  or  range  is 
given  in  Figures  10  and  11.  These  levels  were  obtained  by  coherently 
processing  for  about  0.2  s  and  then  averaging  all  samples  for  one  minute. 
This  procedure  provides  a  sounding  with  a  good  signal-to-noise  ratio. 

Figure  10  was  made  using  the  in-phase  antenna  configuration;  there  is 
no  obvious  antenna  lobe  structure  in  the  in-phase  return;  meteor-caused 
returns  can  be  seen  out  to  200  nmi;  a  small  hump  between  250  and  650  nmi 
is  identified  as  caused  by  sporadic  E;  the  major  return  between  750  and 
1600  nmi  appears  to  be  from  several  one-hop  paths  that  have  washed-out 
antenna  lobe  structure.  The  amplitude-versus-range  pattern  shown  in 
Figure  11  is  that  from  the  anti-phase  antenna  element  feed.  In  it,  a 
peak  at  900  nmi,  a  valley  at  about  1000,  and  another  peak  at  1220  are 
judged  to  be  due  to  the  antenna  pattern.  Using  Figure  8  one  can  calculate: 
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FIGURE  10  RADAR  RETURN,  IN-PHASE  ANTENNA  PATTERN 


i 


m 


p 


<>••••••#•••••  ti*  •  •  *  •  < 


'M 
«r  «t  n 
'  • 
°  r  ■ 
74*0 


5«?.*«SSJ!*5SS22$1} 

M»'<ii4ri»Mi«a««  • 

•  •  .  •  *  •  «  s?  n  ©  •  »  «  •  • 

•  ■'  «••  M  •<!»«»•#« 

•  •  <•  t\  n  «  t  D  <v  <•  «  •  <1  f  • 

•  •«£€.sa»©  «•••«*••  « 

•  •  •  •  t*o<ni»#««  +  «•  «  o 

■  ••OWONMWVi  ©  •  «  • 

•  *  *  4*  •  «>  Jfi  <•  <4  •  •  MS  *  O 
•••••■#»*••*•»*•♦• 

•  •■•««  *IM>  ]«<»••  * 

•  •••••  -v'  • 

••••••••  *>  -«eR»'» 


•••••••  %•  V  o 

•  T  * 

•  •••  i»  ••IUMm'  V  RU 

•  O  •  »  ••  •  •  II  Rij  O  V'  5  « 

•  ••'/•IMMSi-Stl  »v»< 

•  «*  »  •  •••••••It*  • 

•  ••••••••••MS 

•  ••  ••••••••«  <o 


•••••••••••••••••• 54 o 


••••••••••••••••••••a  «o 

••••••••••••••••••«•. « 

<>•••••••••••••••••••00 

•  •••••••••••••••••••«fj 

•  •••••••••••••••••••  )H 

•••••••••••••••••••  • 

•  ••'*•••••■••■••••••  •  * 

••••••••••••••••••••a#  M 

*«••••••••••  •■•  •••••••»  % 

•  •  n  ••••••  •  •  **  o  ••••>  o  4 


•••••••••••^•••••••a  to 

••••••••^^••••••••••♦o 

•  ••••••••••••••••••a  •  w* 

•  ■«*©•••••••••»••••••••  • . 

.••••••••••••••••••••• <4— 

'»«Q««vT  €>••••••  •  •••••••€»•  •  * 

•  •••••  »*••••••••••••«“••••  • 

_••••••<»•••••••••••••••••  • 

•  T/  •  ••••<»  •••••••••••••»o««  *c 

•  •«»•••••  <>••••••••••••••••••5 

•  ••••••••  '»  •••••••••••••••••  •  • 

•  ••••©•••(•••••••••••••••••tp 

_ •  •••••••'•••••••••••••••••••»*» 

••sisj8ss5sss:s::ss:ss::s:sss::0 

••••••••••••••••••••••TO 

•  •••••••••••••••••••••  •  a 

••••••«••••••••••••* 

•••••••••••••••* 

•••••••••••••••  • 

•••••••••«••«••• • 

•••••••••••••••saw • 

•••••••••••••••••••^o 

••••••••••••••••• •“ o 


MIR  M  ••••••••••••••••••••  • 


•  •  •  •  •  •  a  Q 

•  •  •  •  •  •  “  © 


•  ••••••••••••••••••••••••«»•• 

••••••••••••• ••••••••••••••••^ 

•  •••••••*»•••••••••»•••«•••••«.  • 

•  <«••••••••••••••••••••••••••  • 

•  •••»•••••••••••••*>•••••••••  • 

•  ••••  •••••••••••••••••••••  •  a  O 

•  •••••••••••••••••■••••••••a  •  O 

*«•••••••••••••••••••••••••• • • 

•  ••••••••••••••••••••••••a  to 

••••••••••••••••••••••••••a • M 

8  ••••••••••••••••©••••••©•••  •  • 

••••••••••••••••••••••••••a • 

_.••«•••••••••<»••'*«  •»•«••••••••••••  • 

•  ••••••©©••©.••••••••••••••••••••••••» 

:•••••.**  •••••••••••••••©••••••••••mo*  • 

•••••••••••••«•••••••••••••«••••••• ao 

•  •••»•*••••••«»••••••••••••••••••••••••  •  © 

•  ••«•***•••••••••••••••••••••••••••  •  • 

•  •(^••••••••■••••••••••••••••••••O 

•  •••••••••••••••••••••••••a  »  • 

«••••••••••••••»•••••••••• a 

•••••••••••••••••••••••••a • 

•  •©•©••••••••••••••••at)  r 


•assssssssssssssfs 


^••••••••••••••WOM  I 


•  •<»•••••••»••••••••••• • 

•••••••••••••••••••••••••a 

•  •••••••••»  '•••••••••••To 


8  8  8  8 

?  i  2  s 

♦  •  t  • 

<iAjM».aou  mx  o« 


<Z»JH>-3QU  « X  Ofl 


FIGURE  11  RADAR  RETURN,  ANTIPHASE  ANTENNA  PATTERN 


Virtual  Range  (nmi)  Radiation  Angle  (deg) 


900 

22 

1010 

16 

1220 

12 

These  values  are  approximate  and  incomplete  since  they  cover  but  part 
of  the  F-region. 

A  detailed  Doppler  analysis  can  aid  in  path  identification. 
Figure  12  shows  such  analysis  as  6-dB  spaced  amplitude  contours.  The 
upper  plot  is  for  the  in-phase  and  the  lower  for  antiphase  antenna  feed; 
the  ordinate  is  in  Hertz  and  the  abscissa,  R',  in  hundreds  of  nmi.  The 
Doppler  from  sea  returns  at  this  frequency  is  at  about  ±0.4  Hz  and  the 
peak  pairs  with  0.8-Hz  spacing  are  evident.  Land  returns  are  also  seen 
at  a  zero  Doppler.  Low-level  returns  via  sporadic  E  are  identifiable 
between  200  and  600  nmi.  The  in-phase  plot  shows  evidence  of  more  than 
one  path  with  a  multiplicity  of  peaks.  The  antiphase  plot  tends  to  sup¬ 
port  the  prior  tabulation;  further,  the  minimum  in  the  sea  return  at 
R'  =  1180  nmi  can  be  taken  as  identifying  the  widest  part  of  Newfound¬ 
land  at  R  =  1200  nmi.  To  better  unravel  the  in-phase  plot  it  was  re¬ 
analyzed  with  an  expansion  of  the  Doppler  scale,  and  this  is  shown  in 
Figure  13.  The  points  marked  X  are  the  Land-range  estimates  of  equi¬ 
distance  from  the  strait  between  Nova  Scotia  and  Newfoundland.  The  center 
of  the  strait  is  at  R  =  900  nmi  and  for  R'  the  estimates  are  945  and  985 
nmi.  Sporadic  E  is  difficult  to  identify  within  the  larger  returns,  but 
since  it  is  evident  at  the  shorter  ranges  it  is  probably  a  contributor  to 
the  in-phase  returns. 
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FIGURE  13  DOPPLER  ANALYSIS  EXPANDED 


In  summation,  the  transmission-path  description  should  fit: 


R(nmi) 

R/1(nmi) 

R/2(nmi) 

Observed 

900 

945 

985 

Predicted 

900 

Not  seen 

9S2 

Observed 

1120 

1180 

Predicted 

1120 

1184 
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in  addition  to  the  prior  table.  Of  course  the  entire  shape  of-  the  earth 
returns  predicted  from  the  deduced  ionosphere  should  fit  the  observa¬ 
tions.  In  this  particular  case  the  use  of  a  constant-scattering  coeffi¬ 
cient  a0  =  “17  dB  can  cause  some  amplitude  error.  All  of  the  water  illumi¬ 
nated  may  not  have  possessed  this  large  a  scattering  coefficient  and  the 
land  very  probably  did  not. 

It  is  hoped  that  this  exercise  has  shown  the  utility  of  using 
earth  features — here,  the  Doppler  distinctiveness  of  land  and  sea  echoes 
for  path  checkpoints.  This  work  was  done  with  an  azimuthal  beamwidth  of 
10P  ,  one  beam  position,  and  a  receiver  bandwidth  of  1.5  kHz.  Employing 
wider  bandwidth  could  considerably  increase  the  accuracy  of  check-point 
location. 

The  model  was  used  to  generate  median  predictions  for  the 
month  of  April.  These  were  compared  to  the  observed  values  taken  on 
April  30,  1971,  and  revisions  were  made  to  the  median  predicted  ionogram 
to  improve  the  comparison  with  the  measured  backscatter.  Tables  1  and  2, 
and  Figures  14,  15,  and  16  show  the  results  of  running  the  revised  iono¬ 
sphere.  Figure  15,  the  in-phase  return,  shows  the  same  shape  as  measured 
return,  Figure  10.  The  mode  structure,  avowing  combinations  of  E,  Es, 

FI,  and  F2  modes,  eliminates  the  smooth  antenna  pattern  from  the  back¬ 
scatter.  Figure  16,  the  antiphase  return,  however,  shows  the  effects  of 
the  antenna  pattern  because  the  return  is  now  dominated  by  the  F2-layer, 
The  final  check  on  the  analysis  is  to  compare  the  derived  ionogram, 

Figure  17,  with  soundings  taken  near  the  midpath  with  an  ionosonde. 

Figure  17  shows  thi3  check  against  a  vertical-sounder  observa¬ 
tion  near  Boston  that  is  near  the  actual  area  of  reflection.  The 
"adjusted"  ionogram  from  the  one-frequency  backscatter  observations  is 
very  close  to  that  measured  except  for  the  F1-F2  cusp  region.  This  does 
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FIGURE  14  ADJUSTED  IONOG RAM 
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FIGURE  15  ADJUSTED  BACKSCATTER-AMPLITUDES  IN-PHASE  ANTENNA  PATTERN 
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FIGURE  16  ADJUSTED  BACKSCATTER-AMPLITUDES  ANTIPHASE  ANTENNA  PATTERN 
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ACTUAL  OBSERVATION 


not  affect  the  predicted  backscatter  to  any  extent  because  these  rays 
are  highly  absorbed  when  the  deviative  absorption  is  used.  The  remote 
ionosphere  has  been  adequately  sensed  using  one-frequency  backscatter 
data  and  two  antenna  beam  positions. 

The  model  described  in  this  paper  is  a  versatile  tool  that  can  be 
used  in  the  description  of  radar  and  communication  systems  that  use  the 
ionosphere  as  a  transmitting  medium.  The  complexity  of  the  ionospheric 
description  is  intermediate  between  the  simple  four- layer  models5’6  and 
full  ray-trace  programs.  It  is  constantly  being  used  to  predict  aixd 
analyze  observed  data.  Immediate  goals  include  the  construction  of 
methods  to  improve  the  merging  of  the  individual  E,  FI,  and  F2-layers 
(e.g.,  to  determine  when  a  cusp  is  really  present),  and  to  improve  the 
deviative  absorption  factors. 


Appendix 

NUMERICAL  MAPS  OF  fQFl 

Reference  9  describes  the  analysis  of  fcFl  values  for  the 
period  January  1954  through  December  1966.  One  result  of  that  analysis 
is  a  quadradic  fit  of  fQFl  to  the  cosine  of  the  sun's  zenith  angle  and 
a  linear  fit  to  sunspot  number  for  each  month,  defined  as  follows: 

fQFl  =  (A  +  U  ♦  R)  +  (B  +  V  •  R)  cos  (z)  +  (C  +  W  •  R)  cos2  (z)  (1) 

where 

R  =  Sunspot  number 
z  =  Sun' s  zenith  angle 

A,  U,  B,  ¥,  C,  and  W  =  Coefficients  derived  by  a  least  squares 

fit  to  the  data. 


Equation  (1)  is  to  be  used  as  long  as  Z  ^  Z_ov  where  Z  is 

u*ax  max 

given  by  a  linear  fit  to  the  cosine  of  the  Rawer  modified  magnetic  dip 
angle  and  the  sunspot  number  for  each  month,  as  follows: 


Zmax  =  CD  +  X  •  R)  +  (E  +  Y  •  R)  cos  (x) 

where 

R  =  Sunspot  number 

x  =  Rawer  modified  magnetic  dip  angle 

D,  X,  E  and  Y  =  Coefficients  derived  by  a  least-squares  fit 
to  the  data. 
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